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ABSTRACT 

SONYC - Suhstellar Objects in Nearby Young Clusters - is a program to investigate the frequency 
and properties of young substellar objects with masses down to a few times that of Jupiter. Here we 
present a census of very low mass objects in the ~ 1 Myr old cluster NGC1333. We analyze near- 
infrared spectra taken with FMOS/Subaru for 100 candidates from our deep, wide- field survey and 
find 10 new likely brown dwarfs with spectral types of M6 or later. Among them, there are three 
with >M9 and one with early L spectral type, corresponding to masses of 0.006 to < 0.02 M©, so far 
the lowest mass objects identified in this cluster. The combination of survey depth, spatial coverage, 
and extensive spectroscopic follow-up makes NGC1333 one of the most comprehensively surveyed 
clusters for substellar objects. In total, there are now 51 objects with spectral type M5 or later 
and/or effective temperature of 3200K or cooler identified in NGC1333; 30-40 of them arc likely to be 
substellar. NGC1333 harbours about half as many brown dwarfs as stars, which is significantly more 
than in other well-studied star forming regions, thus raising the possibility of environmental differences 
in the formation of substellar objects. The brown dwarfs in NGC1333 are spatially strongly clustered 
within a radius of 1 pc, mirroring the distribution of the stars. The disk fraction in the substellar 
regime is < 66%, lower than for the total population (83%) but comparable to the brown dwarf disk 
fraction in other 2-3 Myr old regions. 

Subject headings: stars: circumstellar matter, formation, low-mass, brown dwarfs - planetary systems 
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1. INTRODUCTION 

Brown dwarfs are objects with masses too low to sus- 
tain stable hydrogen burning (M < 0.08 Mq) and as such 
interme diate in mass between low -mass stars and giant 
planets (jOppenheimer et al.ll2000l) . The substellar mass 
regime is crucial to test how the physics of the forma- 
tion and early evolution of stars depends on object mass, 
thus may help address some of the most relevant issues 
in this field. One example is the origin of the initial mass 
function (IMF) and the relative importance of dynamical 
interactions, fra gmentation, and acc retion in setting the 
mass of objects (|Bonnell et al.|[2007l) . 

Surveys in star forming regions indicate that the mass 
regime of free-floating brown dwarfs extends down to 
mas ses below the Deuterium burning limit at 0.015 
(e.g. IZapatero Osorio et al.l[2000l : E ucas fc R,ochel[2000h . 
i.e. it is overlapping with the domain of massive planets. 
The currently available surveys, however, are not com- 
plete in the substellar regime. Only small regions have 
been observed with sufficient depth to detect the lowest- 
mass brown dwarfs. Moreover, in many cases the brown 
dwarf candidates are selected based oi i their mid-infrare d 
excess and the presence of disks fe.g. lAUers eralll2006f ). 
which introduces an obvious bias. In other cases, the 
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presence of methane absorpt ion is used to identify ob- 
jects fe.g. lBurgess et alll2OO90 . which only finds T dwarfs 
in a limited temperature regime. 

In the SONYC project (short for: Substellar Objects 
in Nearby Young Clusters) we aim for a more complete 
census of brown dwarfs in star forming regions. For 
a number of regions, we have carried out deep photo- 
metric surveys in the optical and near-infrared, to fa- 
cilitate a primary candidate selection based on photo- 
spheric colours. This is complemented by Spitzer data 
to identify additional objects with disks. We have pub- 
lished the photometric survey as well as fol low-up spec- 
troscopy for thre e regions so f ar: N GC1333 (jScholz et al.l 
2009ai) , p Oph dg eerseLaU 1201 1[ ). and Chamaeleon-I 
(jMuzic et aLll201lD . A fourth paper with additional spec- 
troscopy in pOph is submitted (Muzic et al.). Based on 
these results, the largest population of brown dwarfs is 
found in NGC1333, a very young (~ 1 Myr old) cluster 
in the Perseus 0B 2 association at a dista n ce of ^ 300 pc 
(iLad a et al.l Il996t Ide Zeeuw et~all Il999l : iBelikov et all 
2002) . Fig. [1] shows a deep optical image of the clus- 
ter NGC1333 with some of the relevant features marked. 

Here we present new follow-up spectroscopy in 
NGC1333 for a large sample of additional candidates 
from our photometric survey (Sect. [2]) and identify 10 
new, previously unknown very low mass members (Sect. 
[3]). Combining these with the known members yields 51 
objects with spectral type >M5 in this cluster. In Sect. 
|4]we analyse the brown dwarf census for NGC1333, in- 
cluding the mass function, the spatial distribution and 
the disk properties. The conclusions are given in Sect. 
O Throughout this paper, we make use of a large number 
of different samples for objects in the area of NGC1333. 
In Table [1] we provide an overview of the most important 
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TABLE 1 

Overview of the various samples in NGC1333 used in this 

PAPER 



Sample No. 



Objects observed with FMOS (Sect. \^ 100 

excluded as very low mass sources 63 
confirmed as young very low mass sources (YVLM) 26 

newly identified 10 

Objects with spectral type >M5 (Sect. |4J| 51 

identified in this paper 10 

identified in Scholz et al. (2009a) 20 

with estimated masses < 0.08 Mq 30-40 

with Spitzer counterpart (Sect. 14.31 1 41 

with mid-infrared excess at 3-8 fim 27 

Candidates selected from the iz diagram (Sect. I4.1|l 196 

with spectroscopy from MOIRCS or FMOS 98 

confirmed by our spectra 27 

confirmed by other groups 8 

Class I and II sources in NGC1333'' (Sect. 1431 13^ 
with 2MASS detection 94 
with estimated masses 0.02 < M < 0.08 M0 29 
corrected for disk fraction 35 

IGutermuth et al.l t2008i) 




3:30:00.0 30.0 29:00.0 28:30.0 

Right ascension 



Fig. 1. — Subaru/Suprime-Cam i-band image for our target 
region NGC1333 in the Perseus star forming region. Marked are the 
cluster radius with a large dashed circle, the two well-known objects 
BD+30 549 (north) and HH 7-11 (south) with smaU ellipses, and 
the population of very low mass members (see this paper, Tables 
[2] and [4]| with small circles. The reflection nebula NGC1333 is 
slightly north of the image center. Coordinates are J2 000. For 
more information on this image, see IScholz et all ll2009al) . 

samples and link them to the corresponding sections of 
the paper. 

2. OBSERVATIONS AND DATA REDUCTION 

Our spectra were obtained with the Fiber Multi Ob- 
.i ects Spectrograph (FMOS) at the Subaru Telescope 
(jKimura et al.ll20ia) on the night of 2010 November 27. 
Four-hundred fibers, each of 1.2" diameter are configured 
by the fiber positioner system of FMOS in the 30' diam- 
eter field of view, with an accuracy of 0.2" rms. The 
patrol radius of each spine is 87" , while the minimum 
spacing between two neighboring spines is 12" . 

The spectra are extracted by the two spectrographs 



(IRSl and IRS2). Our data were obtained in shared- 
risk mode, using only one of the spectrographs (IRSl, 
200 fibers). IRSl is equipped with a 2048 x 2048 pixel 
HAWAII-II HgCdTe detector, and a mask mirror for 
OH airglow supression. With the low-resolution mode 
(i? ~ 600) the spectrograph yields a coverage from 0.9 
to 1.8 fim (J- and H-bands). We obtained 10 exposures 
with 15 min on-source time each. 

The observations were carried out in the Normal Beam 
Switching (NBS) mode, i.e. the same amount of time was 
spent to observe the sky, which is achieved by ofFseting 
the telescope by 10-15". The seeing during the science 
observations was stable at ~ 1" . Since no stars suitable 
for telluric correction are found within our science field- 
of-view, we observed several standard star fields, covering 
the range of airmasses at which the science target was 
observed (airmass between 1-2). The observed standard 
stars have spectral types in the range F4 - G5. 

Data reduction was carried out using the data reduc- 
tion package supplied by the Subaru staff. The package 
consists of IRAF tasks and C programs using the CFIT- 
SIO library. The data reduction package contains the 
tasks for standard reduction of NIR spectra, perform- 
ing sky subtraction, bad-pixel and flat-field correction, 
wavelength calibration, flux calibration and telluric cor- 
rection. In this last step, each 15-minute exposure was 
calibrated using a standard star at the appropriate air- 
mass. Finally, the ten individual exposures were aver- 
aged. For the analysis (Sect. [3]) the spectra were binned 
to a uniform wavelength interval of 5 A and smoothed 
with a small-scale median filter. For the reduced spec- 
tra, the signal-to-noise ratio in the H-band ranges from 
10 to 70. 

In total, we covered 100 targets, from which 71 are 
selected from our IZ candidate catalogue (jScholz et al.l 
[200 9a) . 10 additional targets have been selected by 
combining our JK-catalogue with the 'HREL' catalogue 
from the Spitzer 'C ores to Disk' (C2D) Legacy program 
(|Evans et al.ll2009[ ). All 10 have colour -excess in IRAC 
bands and thus should have a disk (see Sect. 14. 4p . 
19 fibres were placed o n kno wn M-type members for 
reference (fr om Wilkin g et al.|[20 04: Grei ssl et al.l 120071: 
IScholz et al.l[2009a; Winston et al..,2009. ). 

To test for possible effects of imperfect calibration, we 
compared the H -band spectra for six objects observed 
with MOIRCS (IScholz et al.l 12009a') and with FMOS 
(this paper). For three of them there is excellent agree- 
ment (SONYC-NGC1333-1, 5, 8), while for the others 
there are slight differences in the spectral slope. Using 
the method outlined in Sect. 13.21 we measured the spec- 
tral types for the four out of six MOIRCS spectra which 
cover the entire H-band. All four give types that are 
later than those derived from FMOS, by 0.4, 1.3, 0.5, 0.7 
subtypes, i.e. the differences are larger than our internal 
accuracy of 0.4 subtypes. 

These discrepancies may indicate residual problems 
with the telluric calibration in the MOIRCS and/or the 
FMOS data. These problems could be induced by the 
use of multi-object facilities: Since we stay on the target 
field for long integration times and the fields themselves 
do not cover adequate telluric standards, there is a sig- 
nificant time and position offset between science targets 
and standards, i.e. the depth of the telluric bands could 
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potentially be quite different between science and stan- 
dard fields. 

Stable conditions, as they were present for the FMOS 
observations, should mitigate this effect. The FMOS 
data also have the wider wavelength coverage, which fa- 
cilitates the telluric correction and the spectral analysis. 
Therefore we put more trust in the quantities derived 
from FMOS spectra. 

3. SPECTRAL ANALYSIS 

3.1. Selection of young very low mass objects 

In total we obtained spectra for 100 objects with 
FMOS. Based on the broadband spectral shape in the 
near-infrared, young very low mass sources can be reli- 
ably separated from more massive stars. Objects with 
very low masses and thus effective temperature below 
~ 3500 K or spectral types of ~M3 or later show broad 
absorption features of H2O, which distinguishes them 
clearly from the smooth near-infrared slopes o f more 
massive and hotter stars (jCushing et al.l 120051 ). The 
depth of these absorption troughs is a strong function 
of effective temperature. 

The most important spectral feature for our purposes 
is the H-band 'peak', formed by the two H2O absorption 
bands at 1.3-1.5 and 1.75-2.05 /im. The shape of this 
feature is gravity sensitive; while it appears round with 
a flat top in evolved field objects, it is triangular with a 
pronounced peak at 1. 67 /xw in young objects with ages 
of < 100 Myr ( KirkpatncketjD 120061: iBrandeker et al.l 
I2006t iBihainet al.ll2010l ). In addition. H2O absorption 
causes a sharp edge at 1.35 fim and another 'peak' in the 
K-band. 

We use these features to identify young very low mass 
sources in the FMOS sample. We are looking for objects 
showing structure over the full spectral regime, as op- 
posed to a smooth slope. In particular, we select objects 
with a) a pointy peak in the H-band and b) an absorp- 
tion edge at 1.35 fj,m. Out of 100 FMOS spectra, 26 show 
these characteristics and are called YVLM sample (short 
for 'young very low mass') in the following. For 11 the 
spectra are too noisy to make a decision, and the remain- 
ing 63 do not show these features. These 63 objects for 
which we can exclude that they are very low mass sources 
are listed in Appendix El 

From the 19 literature sources, 16 are re-identified and 
are part of the YVLM sample. The other 3 have spectra 
that are too noisy to identify the features. The remain- 
ing 10 YVLM objects are newly confirmed very low mass 
members of NGC1333 and are listed in Table [5J We use 
the nomenclature SONYC-NGC1333-X for these objects, 
where X is a running number. S ince we have listed ob- 
jects 1-28 in lScholz et al.l (l2009aD. we continue here with 
no. 29. Note that the list in IScholz et al.l ()2009aD con- 
tains some previously confirmed members. The spectra 
for the 10 new objects are shown in Fig. [2j 7 of the 
new objects come from our IZ catalogue, the remaining 
3 from the JK plus Spitzer list (SONYC-NGC1333-31, 
32, 33). 

Many of our spectra show emission in Paschen /? at 
1.28 fim. If the emission originates from the source it- 
self, it would be a clear indication for ongoing accretion 
(jNatta et al.ll2004f ). However, this is difhcuh to tell with 
fibre spectroscopy; our spectra may be affected by emis- 
sion from the cloud material in NGC1333. Still, it is 
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Fig. 2. — FMOS spectra for the 10 new brown dwarfs listed 
in Table [2] The fluxes are on a relative scale with arbitrary off- 
sets. The observed spectrum is shown in black, the best fitting 
DUSTY model slightly offset, with thin, dashed lines. We masked 
out the Pa/3 lines for clarity. The identification numbers SONYC- 
NGC1333-X are abbreviated with S-X. 

worth pointing out that the fraction of objects with Pa /3 
emission is 50% in the YVLM sample and 32% among the 
remaining objects. This indicates that the other sources 
might contain a number of young stars in NGC1333 with 
temperatures > 3500 K. 

In Fig. |3] we compare the FMOS spectra for 
two of our sources, both originally identified in 
IScholz et al.l (l2009aD. with spect ra for young brown 
dwarfs from iMuench et al.l ()2007l, left panel) aiid old 
dwarfs with similar spectral t ypes (jBurgasser et al.|[2008l : 
iBurgasser fc McElwainI 12006'. middle panel). At late M 
spectral types, the more rounded peak of the old objects 
becomes visible (compare the young KPN04 with the 
old LHS2924). However, for earlier types the difference 
between young and old is not significant. The plot also 
illustrates that sufficient signal-to-noise is required to dis- 
tinguish between the two types of objects. Our spectra 
often do not have the necessary quality to make this de- 
cision. From these arguments it follows that contami- 
nation by late M-dwarfs in the field cannot completely 
be excluded. Given the compact nature of the NGC1333 
cluster (see S ect. | 4.3p and the lo w space densities of 
such objects (jCaballero et al.ll2008l) . the contamination 
is considered to be negligible (< 1). 

3.2. Spectral types 
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TABLE 2 

New very low mass members in NGC1333. The IDs SONYC-NGC1333-X are abbreviated with S-X. 



ID 


«(J2000) 


5(J2000) 


i' (mag)'' 


z' (mag)^ 


J (mag)'' 


K (mag)'' 




Av^ 


SpT= 




Other namcs^ 


S-29 


03 28 28.40 


+31 16 27.3 


17.637 


16.675 


14.624 


13.624 








M6.9 


3150 




S-30 


03 28 31.08 


+31 17 04.1 


23.650 


21.235 


16.823 


14.079 


9 


9 


M9.3 


2700 




S-31 


03 29 44.15 


+31 19 47.9 






17.386 


15.290 


6 


5 


~M9 


2300 


Sp 132 


S-32 


03 29 03.21 


+31 25 45.2 






15.798 


13.832 


4 


4 


M7.1 


3200 


MB089, Sp 79 


S-33 


03 29 03.95 


+31 23 30.8 






17.143 


14.932 


7 


5 


M8.3 


2500 


MB0116, Sp 83 


S-34 


03 29 06.94 


+31 29 57.1 


20.839 


19.225 


16.541 


14.774 


4 


4 


M7.0 


2950 


Sp 90 


S-35 


03 29 10.18 


+31 27 16.0 


18.988 


17.711 


15.547 


14.127 


2 


2 


M7.4 


3050 


MB094, Sp 96 


S-36 


03 29 25.84 


+31 16 41.8 


23.392 


21.432 


18.53 


17.07 


2 





~L3 


2250 




S-37 


03 29 30.54 


+31 27 28.0 


17.356 


16.344 


13.808 


12.598 


1 


2 


M6.9 


3200 


MB068 


S-38 


03 29 34.42 


+31 15 25.2 


21.091 


19.654 


17.33 


16.15 


1 


1 


M7.7 


2850 





^ J- and K-band photometry from 2MASS or, for SONYC-NGC1333-36 and 38, from lScholz et al.| (|2009aD 

Calculated from the J- and K-band magnitudes using Equ. 1 
^ Corrected after spectral fitting, sec Sect. 13.31 
^ Estimated using the HPI index as defined in Sect. 13. 2| 
^ Estimated by comparing the spectra to models, see Sect. 13.31 

^ identifiers are from lWilking et al.l II2004I . MBO) and the Spitzor survey bv lGutermuth et al.1 | |2008| . Sp) 



Spectral classification for cool dwarfs in the near- 
infrared relies mostly on the broad H2O absorption fea- 
tures mentioned above. A number of spectral indices 
has been suggested in the literature. Relevant for us are 
the indices that are calibrated for young brown dwarfs 
and thus take into account the fact that the absorption 
features are sensitive to gra v ity. S u ch indices have bee n 
pro posed by [W ilking et alj ()2004[ ). lAUers et al.l (|2007l ). 
and [Weights et al. (200§). We find, however, that these 
schemes are of limited use for our low-resolution and 
(mostly) low signal-to-noise spectra, mainly because the 
wavelength offset between the two intervals from which 
the index is derived is relatively small, e.g. 0.058 for 
the H2O index (Allcrs ct al. 2003) or 0.103 for the 
WH index (Weights et al.ll2009D . 

For our purposes we require a flrrx ratio that maximises 
the contrast between the numerator and the denomina- 
tor. This is achieved by using the flux ratio between 
1.675-1.685 /xm and 1.495-1.505 ^m. The first interval is 
chosen because it corresponds to the position of the H- 
band peak for young objects with M6-M9 spectral type. 
The second interval marks the lowest flux level on the 
blue side of the H-band peak for such objects. 

This index, dubbed HPI for H-peak index, is illustrated 
in the left and middle panel of Fig. |3]where we show liter- 
ature spectra for young and old brown dwarfs in Taurus 
with spectral types M7 and M9.5 compared with FMOS 
data for two of the brown dwarfs in NGC1333. The El- 
band peak clearly is highly sensitive to the spectral type 
in this regime. We expect that this index increases from 
mid M to late M spectral types. The figure also demon- 
strates the advantages of this index at low signal-to-noise 
ratio. 

To calibrate the HPI, we use literature spectra for 20 
young brown dwarfs with spectral types M7-M9.5, for 
which spectra are p ublicly availabl e in the SpeX Prism 
Spect ral Libraries (|Muench et alJ l2007t iLooper et alJ 
|2007| ). Their spectra have been dereddened using the 
reddening law A\ = (A/1.235 ^m)^^-^^ • Aj and the re- 
lation A J — 0.3088 • Ay- The optical extinction Ay is 
derived from their J ~ K colours: 

= [(J-/s:) - (J-if)o]/0.1844 (1) 

These relations are based on the extinction law from 



TABLE 3 

Objects used to calibrate the H-peak index. 



Name^ 


SpT^ 


J (mag)'^ 


K (mag)'^ 


Ay (mag)'* 


HPid 


ITG2 


M7.25 


11.540 


10.097 


2.40 


1.050 


J0444 


M7.25 


12.195 


10.761 


2.35 


1.124 


CFHT6 


M7.5 


12.646 


11.368 


1.51 


1.135 


KPN02 


M7.5 


13.925 


12.753 


0.93 


1.136 


KPN05 


M7.5 


12.640 


11.536 


0.56 


1.137 


CFGH3 


M7.75 


13.724 


12.367 


1.94 


1.184 


J0441 


M7.75 


13.730 


12.220 


2.77 


1.098 


CFHT4 


M7.0 


12.168 


10.332 


4.53 


1.009 


MH04 


M7.0 


11.653 


10.567 


0.47 


1.174 


KPN07 


M8.25 


14.521 


13.271 


1.37 


1.121 


KPNOl 


M8.5 


15.101 


13.772 


1.78 


1.203 


KPN06 


M8.5 


14.995 


13.689 


1.63 


1.185 


KPN09 


M8.5 


15.497 


14.185 


1.69 


1.155 


LRL405 


M8.0 


15.28 


13.91 


2.01 


1.126 


J0457 


M9.25 


15.771 


14.484 


1.56 


1.240 


KPN04 


M9.5 


14.997 


13.281 


3.91 


1.342 


KPN012 


M9.0 


16.305 


14.927 


2.01 


1.261 


J1207 


M8.0 


12.995 


11.945 


0.81 


1.226 


J1139 


M9.0 


12.686 


11.503 


0.99 


1.336 


J1245 


M9.5 


14.518 


13.369 


0.81 


1.306 



^ see SpeX Prism Spectral Libraries: 
http : //pono . ucsd. edu/~ adam/browndwarf s/spexprism/ 
^ spectral types from iMuench et al.| l |2007D for the first 17 and 
'Loopc r ct al.| d 20071 ) for the remaining 3 

^ pho tometry from 2MASS, except LRL405 dPreibisch et al.l [20031 : 
.Luhman et al.||2003l) 

^ derived as described in the text 



ICardelli et al] (|1989t) for Ry = 4.0, which is used 
throughout this paper. We note that Ry varies 
within star forming regions typically from 3 to 5; the 
adopted value is an a verage from the values measured by 
ICardelli et al.l (|1989D for diffuse and dense regions of the 
interstellar medium. It is al so a reasona ble average for 
our target region NGC1333 (iCernis' 

For the reasons outlined in Scholz ct al] ()2009al) we 
use (J — K)o = 1.0. The resulting uncertainty in Ay is 
~ 1 mag. For example it is possible that we overestimate 
Ay due to the presence of K-band excess from a disk. 
This induces an uncertainty of up to 0.04 in the HPI. 
Note that the literature spectral types for the calibrators 
have been dete rmined in the optica l by comparison with 
templates (e.g. iBriceho et al]|2002[ ). The internal accu- 
racy of these optical types is typically ±0.25 subtypes. 
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14000 15000 16000 17000 18000 14000 15000 16000 17000 IBOOO 14000 15000 16000 17000 18000 

Wavelength (A) Wavelength (A) Wavelength (A) 



SpT = -0.84 + 7.66 HPI ^ 



y X 



Fig. 3. — Dereddened spectra of SONYC objects in NGC1333 (solid lines) compared with different comparison spectra (dashed linos). 
Spectra are offset by a constant value for clarity. The spectral range for the HPI is indicated with dashed lines. The SON YC-NGC1333- 
X ide ntification numbers are abbreviated with S-X in all three p anels. Left panel: Youns brown dwarfs in Tauru s l|Muench et al.l 
120071 1. Middle panel: M dwarf spectral standards in the field lIBurgaaser et al.ll2003 : Burgasser fc McElwainI [200^ . Right panel: 
Young substcllar companions (M8 - HIP78530B tL afreniere eiran'20111 ). L2 - 2M1207B dPatience et al.l 120101) . L4 - 1RXSJ1609-2105B 
jLafrcnicre ct al. 2010). In the right panel, we approximate the spectral slope for two of the faintest objects with a linear fit, the full 
spectra are shown in Fig. |2] 

HPI around the correlation is 0.039. As outlined above, 
this can be attributed to the uncertainty in the extinc- 
tion. The HPI is properly calibrated for types M7-M9.5, 
but may also hold for later spectral types as the H-band 
peak increases in strength in the L-type regime (see be- 
low). 

We note that this correlation does not hold for field 
dwarfs. As seen in the middle panel of Fig. |3] these 
objects have flatter H-band peaks which results in lower 
values for the HPI for a given spectral type. 

Using this index we determined spectral types for all 
26 objects in the YVLM sample, after dereddening using 
the corrected extinction derived in Sect. 13.31 9 of them 
have published spectral types. Five of them (MB069, 
MB054, MB077, ASR 29, ASR83) are based on an index 
defin ed in the K-band (jWilking et al.ll2004l: IcTeissl et alj 
I2007t ). The o ther four have t he Sp itzer IDs 131, 104, 
118, and 46 in I Winston et al.l (|2009 ) and were classified 
in the optical. Excluding MB054, ASR29, Sp 104 and 
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M spectral type 

Fig. 4. — H-peak index (HPI) calculated as flux ratio between 
1.675-1.685 /xm and 1.495-1.505 urn p lotted vs. (optical) spectral 
type for published spectra from IMuen ch ct al. (crosses, 2007) and 
[Looper ct al. (plusses, 2007). The correlation after excluding the 
outlier MH04 at M7 is SpT = -0.84 -f 7.66 ■ HPI (dashed lino). A 
typical errorbar is overplotted. 

The calibrators are listed in Tabled 

In Fig. |4]we plot the HPI for the 20 calibrators. The 
plot shows the expected correlation between index and 
spectral type. The one outlier at spectral type M7 cor- 
responds to the object MH04 and deviates in its HPI 
by 3cr from a linear fit. Nothing particular can be seen 
in its spectrum which would explain this discrepancy. 
Based on the H-band peak, MH04 appears to have a 
later spectral type than indicated in the literature. Ex- 
cluding MH04, we derive a correlation of 

SpT = -0.84 + 7.66 -HPI (2) 

with an Icr scatter of 0.37 subclasses. The Ict scatter in 



Sp 118 which have a published spectral type of M6 or 
earlier, for which the HPI is not properly calibrated, the 
HPI types deviate by -0.4, -0.2, -0.2, -0.1, -hO.6 subtypes 
from the published ones, which provides some reassur- 
ance in the usefulness of the HPI. The 10 new objects 
have spectral types of M6.9 or later, classifying them as 
likely brown dwarfs. The spectral types for the new and 
known sources are listed in Table [2] and |4] respectively. 

In the right panel of Fig. [3] we compare FMOS 
data for the four latest-type SONYC objects with spec- 
tra for young ultracool objects with p ublished spectral 
types : the M8 companion to HIP78530 (jLafreniere et al.l 
mV^ . t he L2 companion to the TW Hya brown dwarf 
2M1207 (jPatien ce et al."2010"). and the L4 companion to 
1RXSJ1609-2105B (Lafreniere et al. 2010), ah three with 
ages of 5-10 Myr. Using the HPI, we get spectral types 
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of M8, L2, and L3.5 for these three comparison objects, 
consistent with the hterature values, which indicates that 
HPI could be useful for spectral typing of young early L- 
dwarfs as well. In this plot we approximate the spectral 
slopes for the two faintest objects in our sample with 
linear fits on either side of the H-band peak, to facili- 
tate the comparison. The plot demonstrates that objects 
SONYC-NGC1333-1, 30, 31 are about or later than M8 
and clearly earlier than L2, in line with our classification. 
The object SONYC-NGC1333-36 compares well with the 
L2 and L4 templates, which makes it the coolest object 
discovered in NGC1333 so far. 

We note that the redward slope of the H-band peak 
appears anomalously steep for SONYC-NGC1333-31 and 
36. This is not introduced by the linear fit or the treat- 
ment of the spectra, and it is not seen in the other FMOS 
spectra for NGC1333 or pOph (Muzic et al., subm.), 
which makes a calibration problem unlikely. The effect 
is difficult to explain, as these two spectra are very noisy, 
but it could in principle be a real feature, especially since 
the data is well-matched by the model spectrum (see 
Sect. 13. 3p . Better quality spectra are needed to ver- 
ify the parameters for these objects. Since the HPI is 
defined for the blueward slope of the H-band peak, this 
does not affect the spectral typing. 

3.3. Model fitting 

For the 26 objects in the YVLM sample the FMOS 
spectra were com pared with AMES-DUSTY models 
(jAllard et all l20Qlf ) for low gravity (log.g = 3.5 or, if 
not available, 4.0) low-mass stars®. Using the extinction 
law described in Sect. 13.21 we calculated a model grid 
for Teff = 1500 to 3900 K in steps of 100 K and Ay ^ I 
to 20 mag in steps of 1 mag. The models were binned to 
5 A, the same binsize as the FMOS data. 

The fitting was done in a semi-interactive way. Since 
extinction and effective temperature cannot be deter- 
mined separately with low-resolution, low signal-to-noise 
spectra, we started by adopting the Ay determined from 
the J — K colour, using Equ. (1). For each of the YVLM 
objects, we calculated the following test quantity (O - ob- 
served flux; M - model flux; N - number of datapoints): 



X 



I 

N 



N 

E 



(o - My 

M 



(3) 



This was done for the series of models using the 'pho- 
tometric' Ay] we selected the one with the minimum 
which is typically between 0.005 and 0.05 (with one ex- 
ception with 0.2). This means that the average deviation 
between observed and model spectrum in a given wave- 
length bin of 5 A is in the same range as the noise in the 
observations. 

Usually a few model spectra (2-4) give indistinguish- 
able x; we adopt the average Toff from these best fitting 
models. A visual inspection of the observed spectra with 
models for a range of temperatures shows that clear dis- 
crepancies are visible for temperatures which are > 200 K 
different from the adopted value, i.e. the uncertainty in 
the adopted values is ±200 K. We note that effective tem- 
peratures are necessarily model dependent; our values 

^ downloa ded from 

[http : //pho enix . ens-lyon.fr/Grlds/AMES-Dusty/SPECTRA/ | 




Spectral type HPI 

Fig. 5. — Comparison between the spectral types and effective 
temperatures for the YVLM sample, as determined in this paper 
(Sect. 13.21 and 13.311 . A typical errorbar is overplotted. The solid 
line is a linear fit to the datapoints excluding the 3 outliers marked 
with squares. We al so show the effec tive temperatur e scal es by 
IMentuch etaD 1(2003. dashed line) and ILuhman et~all l[2003. dot- 
ted line). 

should only be interpreted in the context of the AMES- 
DUSTY models. The best fitting models are plotted as 
thin, dashed lines in Fig. [2] for the 10 newly identified 
objects. 

In 6/10 cases this initial fit is already convincing. In 
two more cases, it can be improved by slightly adjust- 
ing Ay by 1 mag, which is within the uncertainty for 
Ay (see Sect. 13. 2p . The two remaining cases give the 
best fit when Ay is changed by 2 mag compared with 
the initial estimate. In Table [2] we list the photometric 
and adjusted value for Ay and the best fit value for the 
effective temperature. 

For SONYC-NGC1333-36, the coolest object in our 
sample, we find an effective temperature of 2250 K, 
which is significantly higher than the published values 
for the two comparison objects shown in Fig. [3 (right 
panel). For 1RXSJ1609-2105B iLafreniere et air(t2010f ) 
find 1800 lb 20 K based on DUSTY and Drift-Phoenix 
(|Helling et al.l 120081 models. For 2M1207B, two inde- 
pendent groups determined 1600 ib 100 K, a gain ba sed on 
comparison with DU STY spec tra (iPatien cc ct al. 120101 : 
iMohantv et"aI1l2007( ). although ^emer et al.i (|201lT ) sug- 
gest that the actual value might be as low as 1000 K. The 
DUSTY spectra for < 2000 K are clearly not in agree- 
ment with our spectrum for SONYC-NGC1333-36. One 
possible explanation is that model fitting for the com- 
parison objects has been done over the full near-infrared 
range, whereas in our case Toff is mostly fixed by the 
shape of the H-band peak. We will revisit this issue in 
Sect. [311 

In Fig. [5| we compare the derived effective tem- 
peratures with the spectral types determined with the 
HPI in Sect. 13.21 The expected trend ~ cooler tem- 
peratures correspond to later spectral types - is visi- 
ble, with only three exceptions. For two of these ex- 
ceptions (SONYC-NGC1333-17 and -31), the signal-to- 
noise ratio in the spectrum is very poor (~ 10), i.e. 
the uncertainties in spectral type and temperature are 
large. For SONYC-NGC1333-17 there are two published 
spectral types of M8 and M8.7, which would move the 
datapoint closer to the general trend. The third out- 
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lier (SONYC-NGC1333-33) has an unusual dip in the 
spectrum around 1.65 /im, which decreases our spectral 
type estimate by ~ 1 subtype, but does not affect the 
model fitting. Excluding these three outliers, the dat- 
apoints are well-approximated by a linear trend (solid 
Hue): Teff = 4191 - 157 x SpT (where SpT corresponds 
to the M subtype). 

The plot also s hows the effective tem perature scales 
bvlMentuch et all |2008, dashed line) and lLuhman et aTl 
(|2003', dotted line), which are derived using the optical 
portion of the spectrum. The two scales agree fairly well 
with each other, although the Mentuch et al. scale is 
an extrapolation and has not been directly calibrated for 
Toff < 3000 K. The trend seen in our data is consistent 
with these lines. Our datapoints are, however, mostly 
above the two lines, indicating that we systematically 
overestimate the temperatures (by ^ 200 K) or the spec- 
tral types (by 0.5-1 subtypes). 

We explored possible reasons for this discrepancy. One 
option is a systematic error in the extinction. Decreas- 
ing Ay by 1 mag for all objects would shift their spectral 
types to earlier types, but only by ~ 0.3 subtypes, not 
sufficient to explain the effect. Moreover, this would also 
increase the best estimate for Tefi by typically 200 K and 
thus cause larger disagreement between datapoints and 
published effective temperature scales. The inverse is 
true as well - increasing Ay would lead to lower T^ff, as 
required, but also to later spectral types. Thus, system- 
atic changes in Ay do not resolve the problem. Given 
the good agreement between our spectral types and liter- 
ature values (Sect. 13. 2p . we suspect that the offset is most 
likely due to a problem with the effective temperatures 
and could indicate issues with the used model spectra. A 
more extensive comparison of the effective temperature 
scales is beyond the scope of this paper. 

3.4. Spectral flux distributions 

To further test the spectral fitting from Sect. 13.31 we 
compare the full photometric spectral flux distributions 
(SFD) for some selected sources, as far as available, with 
the AMES-DUSTY model spectra. We include our own 
photometry in izJK as well as Spitzer/IRAC d ata from 
the C2D 'HREL' catalogue ()Evans et al.ll2009[ ). In Fig. 
[6] we show 3 examples. 

Object SONYC-NGC1333-5 is well-matched by a 
model spectrum for Tcff of 2800-2900 K for Ay = 4mag, 
in line with the parameters derived from the FMOS spec- 
trum (Table |4]). This object shows colour excess red- 
wards of 5 /im, presumably due to the presence of a disk 
(see Sect. g^). The best match for object SONYC- 
NGC1333-30 is obtained for Teff of 2500-2700 with Ay 
of 9-10 mag, i.e. the object may be slightly cooler than 
listed in Table [21 but still within the uncertainty. For 
the coolest object SONYC-NGC1333-36 a good match 
is found for temperatures between 1900 and 2100 K and 
Ay of 0-2 mag. This is again somewhat cooler than the 
estimate given in Table [2] 

In general, when comparing with the full SFD the re- 
sults are similar to the spectral fits to individual wave- 
length bands in the near-infrared, maybe except for the 
regime below 2500 K where the SFD comparison yields 
lower temperatures by 200 K. The comparison also il- 
lustrates that the ideal dataset for a characterisation of 
young brown dwarfs would be a spectrum covering the 



entire near-infrared domain from 1 to 8 [im, thus includ- 
ing five broadband features. 

4. THE BROWN DWARF POPULATION IN NGC1333 

The newly identified very low mass objects in this pa- 
per add to the substantial number in the NGC1333 region 
that have already been confirmed in the literature. In 
Table m we compile all previously spectroscopically con- 
firmed objects with spectral types of M5 or later or effec- 
tive temperatures of 3200 K o r below, from Scholz et al.l 
(l2009al).IWilking et al.l (l200l . ICreissl et al.l (I2007D . and 
IWinston et al.l ()2009f ). Whenever possible, we also re- 
measure d the HPI sp ectral types for the sources iden- 
tified in IScholz et all (l2009a)', based on their MOIRCS 
spectra. In the Table, we list coordinates, photometry, 
spectral types, effective temperatures, and alternative 
names. Adding the 10 objects discovered in this paper, 
the entire sample comprises 51 objects. 

Not all these objects are brown dwarfs; some are very 
low mass stars. Based on the COND , DUSTY, and 
BCAH isochrones (Baraffe et al.l l2003t IChabrier et"all 
l20M iBaraffe et al1ll998h ''. the hvdro gen burning limit at 
1 Myr would be reached at effective temperatures of 2800- 
2900 K, which is foui id to correspond to an (optical) spec- 
tral t ype of M6-M7 (|Luhman et al.ll2003l : iMentuch et"all 
|2008| ). In Table H we made the cut at M5 to include all 
borderline cases as well. Taking this into account, the 
number of confirmed brown dwarfs in NGC1333 is about 
35 ± 5. This is currently one of the largest and best char- 
acterised populations of substellar objects in a single star 
forming region. In the following we will investigate the 
mass function, spatial distribution, and disk properties 
for this sample. 

4.1. The number of brown dwarfs 

Based on our comprehensive spectroscopy, we can put 
some constraints on the total number of brown dwarfs 
in NGC1333 and the mass limits of the current sur- 
ve ys. For this purpose , we use the iz survey presented 
in lScholz et al.l ()2009al ). In Fig. [Tjwe plot the iz colour- 
magnitude d i agram for the 196 candidates selected in 
Sch olz et al.l (|2009a[) . The confirmed brown dwarfs (Ta- 
bles |2| and [4]) , either by us or other groups, are marked 
with squares; all objects for which we have obtained spec- 
troscopy with crosses. Note that the iz candidates are se- 
lected only with a cut in colour and a cut in PSF shape 
to rule out extended objects. No other selection crite- 
ria have been used, i.e. this sample is as unbiased as 
possible. 

We have useful spectra for 98/196 candidates; out of 
these 98, 24 are confirmed by our spectra. In total the iz 
sample contains 35 confirmed objects with Ay < 12 mag. 
Thus, we have a yield of 24/98 (24%) and would expect to 
find 24 more objects among the candidates for which we 
do not have spectra. Since 11 of them (35 minus 24) have 
already been confirmed by other groups, the expected 
number of additional very low mass objects from this iz 
selection is 13. 

The low-mass end of the diagram deserves additional 
discussion. The faintest confirmed brown dwarfs in Fig. 
Elare SONYC-NGC1333-1, 30, and 36 at i = 23.4 - 23.7. 

^ dow nloaded fr om 

|http : //p erso . en s-lyon.fr/lsabelle . baraffe/] 
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Fig. 6. — Photometric spectral flux distributions for three of the coolest members in NGC1333 in comparison with AMES-DUSTY model 
spectra. The disk excess for SONYC-NGC1333-5 is clearly visible. The SONYC-NGC1333-X identification numbers arc abbreviated with 
S-X in all three panels. Errorbars for the fluxes are overplotted, but they arc too small to be visible in some of the blue bands. 



TABLE 4 

Previously confirmed very low mass members in NGC1333. The IDs SONYC-NGC1333-X are abbreviated with S-X. 



ID 


a(J2000) 


(5(J2000) 


J fmag)^ 


K fmag)"'" 


SpT2 


rp 2 


Other names'^ 


S-1 


03 28 47.66 


-f31 21 54.6 


17.55 


15.24 


1^9.2" 


2600", 2800'' 


MB0139 


S-2 


03 28 54.92 


-1-31 15 29.0 


15.990 


14.219 


1X17.9", M6.5'', MS'*, M8.6f 


2850", 2600*' 


ASR109, Sp 60 


S-3 


03 28 55.24 


-1-31 17 35 4 

1 'J ^ -L 1 O O • Tl 


15.090 


13.433 


M7.9'=, M8.2^ 


2900*' 


AST! 38 




f\Q OQ rr. 


_1_Q1 1 fi 0*^ 1 


1 R 1 7 

-LO. -L 1 




A/TQ fi/ 
ivi y . u 






S-5 


no no n A 

03 28 56.94 


+61 20 48.7 


15.362 


13.815 


M7.6 , M6 , iV16.8 


2850 , 2900 


MB(J91, op 66 


S-6 


03 28 57.11 


4-31 19 12.0 


17.24 


15.34 


M7.3'', MS*", M8.0-'' 


3250", 2700*" 


MB0148, ASR64, Sp 23 


S-7 


03 28 58.42 


4-31 22 56.7 


15.399 


13.685 


M6.5^ M7.1'=, M7.7^ 


2800'' 


MBO80, Sp 72 


S-8 


03 29 03.39 


4-31 18 39.9 


15.833 


14.000 


M8.2'', M8.5'', M7.4'^, M8.4^ 


2850", 2600'' 


MB088, ASR63, Sp 80 


S-9 


03 29 05.54 


4-31 10 14.2 


17.072 


15.667 


M8^ M8.6-f 


2600'' 




S-10 


03 29 05.66 


4-31 20 10.7 


17.113 


15.485 




2500'' 


MB0143, Sp 86 


S-11 


03 29 07.17 


4-31 23 22.9 


17.87 


15.62 


M9.2-'' 


(2600)'' 


MB0141 


S-12 


03 29 09.33 


4-31 21 04.2 


16.416 


13.150 




(2500)'' 


MBO70, Sp 93 


S-13 


03 29 10.79 


4-31 22 30.1 


14.896 


12.928 


M7.5^ M7.4'=, M8.1^ 


3000'' 


MB062 


S-14 


03 29 14.43 


4-31 22 36.2 


14.606 


13.035 


M7^ M6.6'=, M7.7-'' 


2900'' 


MB066 


S-15 


03 29 17.76 


4-31 19 48.1 


14.803 


12.988 


M7.5\ M6.5'=, M7.8-'' 


3000'' 


MB064, ASR80, Sp 112 


S-16 


03 29 28.15 


4-31 16 28.5 


13.054 


12.091 


M8.5", M7.5^ M7.5^ M9.1-'' 


2850", 2600'', 2761"= 


Sp 164 


S-17 


03 29 33.87 


4-31 20 36.2 


16.562 


15.481 


M7.5'', M8^ MS.jf 


2600", 2500'' 


MBO140 


S-18 


03 29 35.71 


4-31 21 08.5 


18.50 


16.94 


M8.2f 


2500'' 


Sp 129 


S-19 


03 29 36.36 


4-31 17 49.8 


17.91 


16.38 




2700'' 




S-21 


03 28 47.34 


4-31 11 29.8 


15.484 


12.702 




3100'' 


ASR117 


ASR15 


03 28 56.94 


4-31 15 50.3 


15.056 


13.461 


M7.4^ Me.O" 






ASR17 


03 28 57.15 


4-31 15 34.5 


15.405 


13.186 


M7.4'=, M6.0'' 




Sp 68 


MB073 


03 28 58.24 


4-31 22 09.3 


16.004 


13.367 


M6.4'= 




Sp 70 


ASR24 


03 29 11.30 


4-31 17 17.5 


13.977 


12.915 


M8.2'=, MS.O"^ 






MB069 


03 29 24.45 


4-31 28 14.9 


14.041 


12.686 


MT.O", M7.4= 


3200" 




ASR29 


03 29 13.61 


4-31 17 43.4 


16.441 


13.028 




3250" 




ASR105 


03 29 04.66 


4-31 16 59.1 


15.550 


12.665 


M6'' 




Sp 84 


ASR8 


03 29 04.06 


4-31 17 07.5 


13.310 


12.313 


M7'* 






MB078 


03 29 00.15 


4-31 21 09.2 


16.466 


13.349 


M5'^ 




Sp 75 


Sp 45 


03 28 43.55 


4-31 17 36.4 


12.219 


10.138 


M5.0'= 


3125"= 


ASR127 


Sp 46 


03 28 44.07 


4-31 20 52.8 


14.245 


12.627 


M7.3'', M7.5'= 


3050", 2829"= 




Sp 49 


03 28 47.82 


4-31 16 55.2 


12.940 


10.909 


M8.0'= 


2710"= 


ASRlll 


Sp 53 


03 28 52.13 


4-31 15 47.1 


13.161 


12.029 


M7.0'= 


2846"= 


ASR45 


Sp 55 


03 28 52.90 


4-31 16 26.4 


13.616 


12.476 


M5.0'= 


3154"= 


ASR46 


Sp 58 


03 28 54.07 


4-31 16 54.3 


13.027 


11.599 


M5.0'= 


3098*= 


ASR42 


Sp 94 


03 29 09.48 


4-31 27 20.9 


14.154 


12.692 


M5.0'= 


3098"= 


MBO60 


Sp 105 


03 29 13.03 


4-31 17 38.3 


15.231 


14.158 




2710= 


ASR28 


Sp 131 


03 29 37.73 


4-31 22 02.4 


13.987 


12.958 


M7.6'', M7.0'= 


2850", 2891*= 


MB065 


Sp 157 


03 29 12.79 


4-31 20 07.7 


14.676 


13.294 


M7.6'', M7.7=, M7.5'= 


2950", 2812*^ 


MB075, ASR83 


Sp 177 


03 29 24.83 


4-31 24 06.2 


14.433 


13.383 


M6.5'', M6.7'=, M6.5'= 


3000", 2957'= 


MB077 


Sp 71 


03 28 58.24 


4-31 22 02.1 


14.912 


12.406 


M6<= 


2990'= 


MB047 



1 photometry from 2MASS (most objects) or SONYC (for SONYC-NGC1333-1, 4, 6, 11, 18, 19) 

spectral types or effective temperatures from the following references: a - this paper, b -IScholz et al.l ll2009al) . c - IWilking et al.l i fSool) . d - IGreissl et al.l 
||2007f ). e - IWinston et al.l ll200g|), f - HP I spect ral types for spectra i n IScholz etaU 1120093) 

^ identifiers are from lAspin et al.l | |1994 ASR). IWilking et all | |20041 MBO), and the Spitzer survey bv lGutermuth et al.l 1 120081 Sp) 
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Fig. 7. — Colour-magnitude diagram for the iz candidates (dots), 
originally identified in Scholz ct al. (2009a). Crosse s are objects 
for whi ch we have obtained spectra in this paper or IScholz et al.l 
Il2009al ). Confirmed very low mass objects from this paper or 
the literature are marked with squares. The horizontal dashed 
line shows the com pleteness limit of the survey estimated in 
IScholz et aLl Il2009al) . 

Comparing their effective temperatures with the 1 Myr 
DUSTY and COND isochrones, it seems hkely that 
they have masses of < 0.02 Mq. If our estimate of 
Tcff = 2000 K for SONYC-NGC1333-36 is correct (Sect. 
13. 4p . the best mass estimate would be in the range of 
0.006 Afg. 

We have taken spectra for 7 fainter objects, but none of 
them is a brown dwarf, which might indicate that we have 
re ached the 'bo t tom of the IMF', as prehminarily stated 
in IScholz etal\ ()2009aD . However we may not be 100% 
complete in this magnitude range. The formal complete- 
ness limits of the iz survey at i = 24.8 mag, determined 
from the peak of t he histogram of the magnitudes (see 
IScholz et all [2009a'). is shown with a dashed line. This 
limit has been derived for a field of view of 34' x 27', but 
most of the cluster members are located in a smaller re- 
gion of 10' X 10' which is partially affected by significant 
background emission from the cloud (see Fig. [T] Thus, 
the completeness limit in the relevant areas might not 
always reach the value shown in Fig. O 

Thus, based on our new da ta we retract the previous 
claim by IScholz et al.l (|2009aD stating a deficit of objects 
with M < 0.02 Mq in NGC1333, for two reasons: a) The 
updated brown dwarf census contains a few of objects 
with masses at or below 0.02 Mq, including one with an 
estimated mass of 0.006 MQ. b) The current survey may 
not be complete at the lowest masses, i.e. we cannot 
exclude the presence of a few more objects with M < 
0.02 Mq. 

The census for M > 0.02 Mq is more robust. From the 
35 confirmed members in the iz diagram we subtract 10 
which are probably slightly above the substellar bound- 
ary (see discussion in Sect. IH). We also subtract the 3 
which are likely below 0.02 M© and add the 13 which 
we are still missing. This gives a total number of ~ 35 
brown dwarfs down to 0.02 Mq and with Ay ^ 12 mag. 

4.2. The star vs. brown dwarf ratio 

As a proxy for the shape of the mass function, pre- 
vious authors have used the ratio of stars to brown 
dwarfs, where these two groups are defined by a range of 



masses. These ratios are more robust ag ainst uncertain- 
ties in the masses than a complete IMF. lAndersen et al.l 
(|2008D use a range of 0.08-1.0 M© for stars and 0.03- 
0.08 M© for brown dwarfs, hereafter called Other 
authors use 0.08-10 M© for stars and 0.0 2-0.08 M© for 
brown dwarfs, hereafter c alled R2 (e.g. Bricefio et al.l 
I2002t iMuench et all l2002t iLuhmaii et al]|2003D . slnce 
the number of high-mass stars is small, the two ratios 
i?i and i?2 should be fairly similar. 

The comparison with NGC1333 is complicated by the 
fact that no comprehensive spectroscopic census is avail- 
able for the stars. The best starting point is p robably 
the Spitzer analysis bv iGutermuth et al.l ()2008l ). They 
find a total of 137 Class I and II members with disk, 
from which 94 are in 2MASS. Objects not detected in 
2MASS are likely embedded sources with high extinction 
Ay > 10 mag and thus not comparable with our brown 
dwarf sample. We calculated absolute J-band magni- 
tudes for this sample using the dereddening described 
in Sect. 13.21 and assuming a distance of 300 pc, which 
gives a range of Mj = — 9 mag. Compar ing with the 
BCAH IMyr isochrone (jBaraffe et al.lll998[ ) the sample 
contains 13 objects with M > 1.0 M©, 52 objects with 
0.08 < M < 1.0 Mq and 29 with 0.02 < M < 0.08 M©. 
These are only objects with disks; correcting for a disk 
fracton of 83% (IGutermuth et al. 2008) shifts the num- 
bers to 16, 63, 35. The latter number is consistent with 
the estimate of brown dwarfs in this cluster given in Sect. 
14.11 Out of 35 brown dwarfs, the number of objects with 
masses above 0.03 M© would be 28. 

Based on these estimates the ratios for NGC1333 be- 
come Ri = 63/28 = 2.3±0.5 and R2 = 79/35 = 2.3±0.5 
(see below for an explanation of the uncertainties) . Our 
value fo r j?i is somewhat la rger than our first estimate 
given in lScholz et al.l (|2009aD of 1.5 ±0.3, mainly because 
we use here the cutoff at 0.03 M© to be consistent with 
lAndersen et~all ()2008[ ). 

The uncertainties for i?i and R2 stated above are la 
confidence intervals and h ave been derived based on the 
prescription provided by ICameronI ()2011[ ). This pre- 
scription is given for population proportions ('success 
counts'). Therefore, we use the Cameron equation to cal- 
culate the confidence intervals aki for the ratio of number 
of stars to the sum of stars and brown dwarfs (kl) and 
crfe2 for the ratio of number of brown dwarfs to the same 
sum (/c2). The confidence intervals for Ri and R2 are 
then derived as follows: 



kl 



k2 - ak2 
kl - (Tfci 

k2 + (Tfc2 



(4) 



(5) 



In Table [5] we compare the ratios for NGC1333 with 
the available literature values for i?i and R2 for other re- 
gions. The same numbers are plotted in Fig. [Sjfor illus- 
tration. To have accurate and consistent confidence in- 
tervals, we re-calculated the errors for all literature values 
as described above. NGC1333 has the lowest ratios mea- 
sured so far in any star forming region, suggesting that 
the number of brown dwarfs in NGC1333 is unusually 
high, which is in line with the conclusion in lScholz et al.l 
([2009a). In particular, the ratios for NGC1333 deviate 
by more than 2a from those in IC348. It should be noted 
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TABLE 5 

Star, to brown dwarf ratios for various star 
forming regions (sbe text) 



Region 








NGC1333 


2.3'> 


(1.8-2.8) 


2.3'' (1.8-2.7) 


ONC 


3.3<= 


(2.8-3.9) 


3.8,4.5,5.0f'S.i> 


UpSco 


3.8<^ 


(3.1-4.5) 




NGC2024 


3.8<= 


(2.6-5.2) 




Chamaeleon 


4.0<= 


(2.3-6.0) 


3.9' (2.9-5.0) 


p-Oph 


5.1" 


(3.8-6.4) 


4.8"= (3.7-6.0) 


Taurus 


6.0<= 


(4.5-7.7) 


7.3^ (5.1-9.6) 


IC348 


8.3<= 


(6.4-10.5) 


8.0^ (6.3-10.0) 



^ In brackets la confidence intervals, see text. Note that 
for tlie i?2 value in the ONC no absolute numbers arc 
available, thus no error estimate is possible. For p-Oph 
we adopted the average numbers from the ranges given by 
Muzic et al. (subm.). 



this paper 



^ [Andersen e t al.| j200gD 

IPawson ct aEI i2onh 
^ Muzic et al. (subm.) 
f [Liihman ct al.' f2003f) 
S|Muoncli ct al. (2002) 
^ [Slcsniclc et al.. 12004 ) 
' ILuhmanI ll2007f) 



I — 

1 - 



o o o 



NGC1333 
ONC 
UpSco 

NGC2024 
, ' Cham 
" , ' RhoOph 



Taurus 



IC348 



star vs. brown dwarf ratio 



Fig. 8. — Star vs. brown dwarf ratios for various star forming 
regions, from Table [5] Ri is plotted with squares, R2 with circles. 
The number on the y-axis is identical with the row in Table |5] 



that the current census for IC348 (iLuhman et all 120031) 
is nearly complete down to 0.03 and covers most of 
the cluster, which makes the difference to NGC1333 even 
more striking. 

This finding has to be substantiated with further sur- 
vey work in diverse regions. In Table [S] we list the statis- 
tical la confidence intervals, purely based on the sample 
sizes. These statistical errors do not take into account ad- 
ditional sources of uncertainty, e.g. unrecognised biases, 
inconsistencies in sample selection or problems with the 
mass estimates, i.e. the actual errors may be larger than 
listed in Table [3 In particular, it is important to note 
that all mass estimates are necessarily mo del- dependent. 
For the value in NGC1333 we use the BCAH isochrones, 
mainly because they cover the brown dwarf regime down 
to the Deuterium burning limit. The problems and un- 
certainties of thes e type of models at v ery young ages are 
well-documented (jBaraffe et al.|[200^ . 



X Class I/II 
a BDs/PMOs 




52.5 52.4 52.3 52.2 52.1 52 

RA (deg) 

Fig. 9. — Spatial distribution of NGC1333 members. Crosses 
are all 137 objects with Spitzer excess by [Gutermutli et al. (2008), 
squares are all confirmed brown dwarfs (Tables [2] and [4]| . Objects 
with spectroscopy for which we can exclude that they are substellar 
members are shown with dots. 

The best way of assessing the overall uncertainties is 
to compare results from independent surveys. As can 
be seen in Table [5l so far the results from independent 
groups agree within the statistical errorbars, with the 
possible exception of the ONC.^ Such an independent 
confirmation is required for NGC1333 as well. 

If confirmed, the unusually low ratio of stars to brown 
dwarfs in NGC1333 could point to regional differences 
in this quantity, possibly indicating environmental dif- 
ferences in the formation of very low mass objects. One 
option to explain this is turbulence, as very low mass 
cores which can potentially collapse to brown dwarfs 
could be a ssembled by the turbulent flow in a molec- 
ular cloud (jPadoan fc Nordlundl[2004D . At first glance 
this could be a realistic possibility for NGC1333, where 
the cloud is strong ly affected by numerous outflows 
(jQuillen et al.l l2005l). although it is not clear if the tur- 
bulence in NGC1333 is mainly driven by these outflows 
(jPadoan et al.l l2009l ). Alternatively, additional brown 
dwarfs could form by gravitatio nal fragmentation o f 
gas falling into the cluster center (jBonnell et al.l 120081 ). 
This latter mechanism would benefit from the fact that 
NGC1333 has a higher stellar density and thus a stronger 
cluste r potential than mos t other nearby star forming re- 
gions (IScholz et al.ll2009aD . 

4.3. Spatial distribution 

In Fig. [9] we show the spatial distribution of the sam- 
ple of very low mass objects listed in Tables [5] and H] 
(squares). For compari son, the positions of the 137 Class 
I and Class II sources (jGutermuth et all 120081 ) are over- 
plotted with crosses. The dots indicate the positions of 
all targets for which we have obtained spectra but which 
are not confirmed as very low mass objects. Addition- 
ally, we show the frequency of objects as a function of 
distance from the cluster center in Fig. [TUl again for the 
same three samples, and in addition for all photometric 
candidates from our iz catalogue (dash-dotted line). 

In the two figures, the spatial distribution of brown 



^ A new paper by lAndersen et al.l Il2011h updates the value of 
for the ONC to Ri = 2.4 ± 0.4, based on an HST survey covering 
a larger area than previous studies. 
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BDs/PMOs 
Class I/II 
no young BDs 
candidates 



0.2 0.3 0.4 

Distance (deg) 

Fig. 10. — Histogram of the distance from the cluster center for 
confirmed brown dwarfs (solidline. Tables [2] and , objects with 
Spitzer excess by Gutcrmu th et aLl || 2008 ) (dashed line), objects 
for which we can exclude that they are substellar members (dotted 
line), and all IZ candidates (dash-dotted line). As cluster center 
we used the average position of the Class I/II sources. 



dwarfs is strongly clustered and indistinguishable from 
the distribution of the total population of Class I/II 
sources in NGC1333. For the two samples, the average a 
and i5 differ only by 0.4' and 0.3', respectively, which is 
< 10% of the cluster radius. Adopting the average posi- 
tion of the Class I/II sources as cluster center, the aver- 
age distance from the center is similar in the two samples, 
5.2' for the very low mass objects and 5.5' for the Class 
I/II sources. The fraction of objects with distance from 
the cluster center of d < 0.1, 0.1 < d < 0.2, d > 0.2 deg 
is 65, 35, 0% for the very low mass objects and 67, 26, 
6% for the Class I/II objects. The scatter in the posi- 
tions is (Ja = 0.071 and as = 0.069 deg for the very low 
mass objects, and ctq, — 0.067 and as = 0.085 deg for the 
Class I/II sources. For all these quantities there are no 
significant differences between the two samples. 

The figures also show that our spectroscopic follow- 
up covers an area that is about 1.5-2 as large (in radius) 
than the cluster itself. We took spectra for 31 candidates 
with distances of > 0.2 deg from the cluster center, but 
none of them turned out to be a brown dwarf. There 
are still 43 candidates from the IZ photometry outside 
0.2 deg (see Sect. I4.ip for which wc do not have spectra, 
but based on our current results, it is unlikely that they 
contain any very low mass cluster members. Thus, our 
wide-field follow-up spectroscopy shows that there is no 
significant population of brown dwarfs at d > 0.2 deg 
from the cluster center, corresponding to ^ 1 pc at a 
distance of 300 pc. 

It has been suggested that gravitational ejection oc- 
curs at an early stage in the evolution of substellar 
objects, either from multiple stellar /substellar systems 
(|Reipurth fc C larke 2001; Umbrcit ct al. 20 0^) or from 
a protoplanetary disk ( Stamatellos &: Whitworthll2009[ ). 
This ejection is thought to remove the objects from their 
accretion reservoir and thus sets their masses. In these 
scenarios one could expect the brown dwarfs to have high 
spatial velocities in random directions. 

An ejection velocity of Ikms"^ would allow the ob- 
ject to travel 1 pc in 1 Myr, i.e. in the case of NGC1333 
this would be sufficient to reach the edge of the cluster. 



However, the gravitational potential of the cluster will 
significantly brake the motion of the brown dwarf. As- 
suming a total cluster mass of 500 Mq (iLada et al.lH"996l) 
homogenuously distributed in a sphere with 1 pc radius, 
a brown dwarf that gets ejected in the cluster center with 
1.5, 2, 3kms"^ would reach a velocity of approximately 
0.5, 1.4, 2.6kms~^ at a distance of Ipc from the cen- 
ter. All objects with ejection velocities of > 2kms~^ 
would have moved to distances significantly larger than 
1 pc over 1 Myr. As shown above, the presence of such 
objects can be excluded f rom our data. 
The scenarios by Umbreit et al.l ()2005f ) and 



IStamatellos fc WhitworthI (|2009( ) predict that a sub- 
stantial fraction of ejected brown dwarfs (more than 
50% in some simulations) exceed this velocity threshold 
of 2kms^^. These models would require some tuning 
to reproduce a spatial distribution as observed in 
NGC1333. However, such simplified scenarios do not 
take into account that dynamical interactions affect 
the total cluster population, not exclusively the b rown 
dwarfs. The cluster formation simulations by iBatd 
()2009() show that the velocity dispersion in a dense 
cluster is not expected to increase in the very low mass 
regime. Although brown dwarfs undergo ejection in 
the simulations, this does not lead to a velocity offset 
in comparison to the stars. NGC1333 seems to be 
consistent with this picture. 

As a side comme nt, we note t hat the parameters in the 
main simulation in iBatd (|2009t ) with gas mass of 500 Mq 
and cloud radius of 0.4 pc are fairly similar to the prop- 
erties of NGC1333, although the simulation produces a 
much higher number of stars and brown dwarfs (total 
stellar mass of 191 Mq vs. ~ 50 Mq in NGC1333). 

4.4. Disks 

In Fig. [TT]we plot the Spitzer/IRAC colour-colour dia- 
gram for the sample listed in Tables [2] and IH again based 
on the C2D-'HREL' catalogue. Out of the sample of 51 
sources with confirmed spectral type M5 or later, 41 have 
photometric errors < 40% in all four IRAC bands and are 
shown in this plot. The figure shows the typical appear- 
ance with two groups, one around the origin, the second 
with significant colour excess in mid-infrared bands due 
to the presence of circum-(sub)-stellar material. In this 
sample of 41 objects, 27 show evidence for a disk, i.e. 
66 ± 8%. All of them have colours comparable to the 
Class II sources identified in lGutermuth et al.l ()2008D . 

The derived disk fraction of 66% is only valid for 
the sample of 41 objects with reliable Spitzer detec- 
tion. In the entire sample of 51 very low mass sources 
in NGC1333 listed in Tables H and El the disk fraction 
is likely to be smaller, because the ten objects which are 
not detected by Spitzer are unlikely to have a disk. Cor- 
recting for this effect, the disk fraction in the full sample 
could be as low as 27/51 or 55%. Therefore, we consider 
the disk fraction of 66% to be an upper limit. 

For compariso n , for the total cluster population 
iGutermuth et al.l ()2008[ ) derive a disk fraction of 83% 
from a Spitzer survey. This number has been derived for 
object s with K < 14 mag. This magnitude limit was cho- 
sen bv lGutermuth et 311(120081 ) because it corresponds to 
M = 0.08 Mq at age of IMyr and Ay = 20 mag. Their 
sample thus includes mostly stars, but also some brown 
dwarfs (as evident from Table |4l which contains a num- 
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0.5 1 1.5 

[5.8]-[8.0] 

Fig. 11. — Spitzer/IRAC colour-colour diagram lor the very low 
mass sources listed in Tables[2]and|4] using the 'HREL' photometry 
from the C2D survey. We only plot objects with photometric errors 

< 40% in all four bands (41 out of 51). The position of the Class 
II objects from the survey bv lGutermuth et al.l I I2008I ) is indicated 
by the dashed box. Objects with reliable detection at 24 fj,m are 
marked with squares 

ber of objects from the lGutermuth et al.l (|2008| ) sample, 
marked with 'Sp'.) The Spitzer sample contains a sub- 
stantial number of objects with Ay > 10 mag, which are 
rare among the currently known brown dwarfs. It is pos- 
sible that some of the heavily embedded brown dwarfs 
with Ay > 10 mag have not been found yet. This could 
explain the discrepancy in the disk fractions. 

Our disk fraction is consistent with the values de- 
rived for ver y low mass members in a Ori, Chamaeleon-I, 
and IC348 (jLuhman et al.l I2008D although all three re- 
gions are thought to be somewhat older (2-3 Myr) than 
NGC1333. 

A more detailed SED analysis was carried out for ob- 
jects with an additional datapoint at 24 /xm. 19 of the 
objects in Fig. |ll]have MIPS fluxes at 24 /Ltm with errors 

< 40%. At this wavelength the images are strongly af- 
fected by the cloud emission and blending. To make sure 
that the fluxes are trustworthy, we checked all objects in 
a 24 /xm image obtained in the Spitzer program #40563 
(PI: K. Wood, AOR 23712512), which is deeper than the 
C2D mosaics. After visual inspection, 3 objects were dis- 
carded; the remaining 16 are point sources at 24 fim and 
are marked with squares in Fig. 1111 

In Fig. [12] we show their SEDs after dereddening (see 
Sect. 13 2p and scaling to the J-band flux (crosses). For 
comparison, the photospheric SED from a model spec- 
trum is overplotted with small dots. To assess the disk 
evolution in the substellar regime, we derive the typical 
SED for NGC1333 and three other star forming regions: 
pOph (IMyr), Taurus (2 Myr), and UpSco (5 Myr). For 
this purpose we selected the objects which are detected in 
all four IRAC bands and at 24 /urn. For p Oph we started 
with the census in Muzic et al. (subm.) and made use 
of the C2 D data. For Taurus we us ed published Spitze r 
data from iGuieu et al.l (l2 p07D and [S cholz et alj (|2006l ). 
For UpSco the data from lScholz et a l. (2007) was used. 
When comparing the SEDs from different regions, one 
has to take into account that the depth of the 24 /im 
observations is not the same; thus the median SED is 
affected by incompleteness at low flux levels. Instead we 
plot the SED for the object that has the 10th highest flux 




Wavelength (/^m) 



Fig. 12. — Spectral energy distributions for 16 very low mass ob- 
jects with disks in NGC1333 (crosses). The SEDs have been dered- 
dened and scaled to the J-band fluxes. We overplot the typical SED 
for NGC1333 (solid line), pOph (dotted line) Taurus (dashed line), 
and UpSco (dash-dotted line). We also sh ow the photospheri c SED 
from the 2800 K DUSTY-AMES model l|Allard et al.ll200ll ) with 
small dots. A typical errorbar for the 24 fi fluxes is overplotted. 

level at 24 iim after converting to XFx and scaling to the 
J-band flux. This represents an estimate for a typical 
SED unaffected by the depth of the Spitzer observations 
and the distance to the cluster. Note that all objects used 
for this exercise are spectroscopically confirmed members 
of the respective clusters. 

For wavelengths < 8 fim the four median SEDs are 
fairly similar. At 8 /im the SEDs in the youngest regions 
(NGC1333, pOph) are slightly enhanced. The biggest 
differences are seen at 24 /xm, particularly when compar- 
ing NGC1333 with UpSco. This is mostly due to the 
fact that NGC1333 harbours a few objects with unusu- 
ally strong excess emissio n, which are missin g in UpSco 
(compare with Fig. 1 in fScholz et al.l [20071 ). indicating 
that the objects in NGG1333 are in an early evolutionary 
st age compared with th e other regions. As demonstrated 
in lScholz et al.l (|2009b[ ) a large spread in 24 fj.m fluxes, as 
seen in NGG1333, can easily be explained by a range of 
flaring angles in the disks. 

5. CONCLUSIONS 

As part of our survey program SONYC, we present 
a census of very low mass objects in the young clus- 
ter NGC1333 based on new follow-up spectroscopy from 
Subaru/FMOS. To derive reliable spectral types from our 
data, we define a new spectral index based on the slope 
of the H-band peak. We find 10 new likely brown dwarfs 
in this cluster, including one with a spectral type '~L3 
and two more with spectral type around or later than 
M9. These objects have estimated masses of 0.006 to 
0.02 Mq, the least massive objects identified thus far in 
this region. This demonstrates that the mass function 
in this cluster extends down to the Deuterium burning 
limit and beyond. By combining the findings from our 
SONYC survey with results published by other groups, 
we compile a sample of 51 objects with spectral types 
of M5 or later in this cluster, more than half of them 
found by SONYC. About 30-40 of them are likely to be 
substellar. The star vs. brown dwarf ratio in NGC1333 
is significantly lower than in other nearby star forming 
regions, possibly indicating environmental differences in 
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the formation of brown dwarfs. We show that the spatial 
distribution of brown dwarfs in NGC1333 closely follows 
the distribution of the stars in the cluster. The disk 
fraction in the brown dwarf sample is < 66%, lower than 
for the stellar members, but comparable to the brown 
dwarf disk fraction in 2-3 Myr old regions. The substel- 
lar members in NGC1333 show a large fraction of highly 
flared disks, evidence for the early evolutionary state of 
the cluster. 
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APPENDIX 

SPECTROSCOPICALLY EXCLUDED OB.JECTS 

In Tables [6] and [7] we provide a full list of objects for which we obtained spectra and which were not classified as 
young very low mass objects based on the shape of their near-infrared spectrum (see Sect. 13. ip . The spectra come 
from the first campaign with MOIRCS (jScholz et al.ll2Q09al ) and from the second run with FMOS (this paper). Most of 
these objects are likely to be either young stellar objects in NGC1333 or background stars with efi'ective temperatures 
above ^ 3500 K or spectral type earlier than ~M3. In Table[n]we alsogive t he J- and K-band pho tometry from 2MASS 
and the identifiers from the photometric surveys by Lada et al. (199l| and lWilking et ahl (|2004l) . if available. Objects 
without listed identifiers are not known to have a counterpart within 1" . 
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TABLE 6 

Objects excluded by spectroscopy in this paper 



a(J2000) <5(J2000) J (mag)l K (mag)l Sel^ Spec^ Idontifior* 



03 


29 


18.71 


+31 


32 


26.4 


16.722 


13.839 


IZ 


F 






03 


29 


52.35 


+31 


28 


13.7 


14.919 


13.829 


IZ 


F 






03 


29 


48.12 


+31 


28 


29.4 


15.550 


13.663 


IZ 


F 






03 


29 


37.80 


+31 


27 


48.4 


17.274 


15.249 


IZ 


F 






03 


29 


34.76 


+31 


29 


08.1 


13.647 


11.532 


IZ 


F 


MBO 43 




03 


29 


11.53 


+31 


30 


05.6 


16.740 


13.786 


IZ 


F 


[LAL96] 241, 242 


03 


28 


55.74 


+31 


30 


58.0 


15.175 


13.390 


IZ 


F 


LAL96 143 


03 


28 


52.66 


+31 


32 


04.3 


16.563 


14.730 


IZ 


F 






03 


28 


46.67 


+31 


31 


35.4 


15.647 


13.933 


IZ 


F 






03 


28 


37.55 


+31 


32 


54.5 


15.107 


14.232 


IZ 


F 






03 


28 


23.84 


+31 


32 


49.3 


14.714 


13.655 


IZ 


F 






03 


28 


44.96 


+31 


31 


09.9 


17.057 


15.211 


IZ 


F 






03 


28 


46.24 


+31 


30 


12.1 


15.089 


14.069 


IZ 


F 


Ft A T Oi^l 


QQ 
OO 


03 


29 


45.53 


+31 


26 


56.5 


14.870 


14.023 


IZ 


F 




03 


30 


05.41 


+31 


25 


13.1 


14.921 


13.680 


IZ 


F 






03 


28 


49.48 


+31 


25 


06.6 


15.957 


13.658 


IZ 


F 


"XTTy/^ on 




03 


29 


10.46 


+31 


23 


34.8 


15.633 


12.762 


JK 


F 


Ft A T Ofil 




03 


28 


43.17 


+31 


26 


06.1 


15.617 


13.384 


IZ 


F 


Ft A T rii^l 


78 


03 


28 


37.75 


+31 


26 


32.8 


16.847 


14.340 


IZ 


F 


[LAL96] 


bU 


03 


27 


56.27 


+31 


27 


00.8 


16.836 


13.967 


IZ 


F 






03 


28 


07.64 


+31 


26 


42.4 


15.974 


13.719 


IZ 


F 






03 


28 


19.51 


+31 


26 


39.5 


15.100 


13.687 


IZ 


F 


Ft A T Qd] 


1 Q 

io 


03 


28 


40.22 


+31 


25 


49.1 


15.637 


12.911 


IZ 


F 






03 


28 


47.64 


+31 


24 


06.2 


14.199 


11.660 


JK 


F 


Ft A T 0(\] 


yo 


03 


28 


55.22 


+31 


25 


22.4 


14.735 


12.550 


IZ 


F 


Ft A T Oi^l 




03 


29 


03.32 


+31 


23 


14.8 


17.254 


14.071 


JK 


F 


Ft A T fltil 


191 


03 


29 


03.13 


+31 


22 


38.2 


13.724 


11.323 


JK 


F 


Ft A T Ofil 
[L;AijyOj 




03 


29 


27.61 


+31 


21 


10.1 


14.836 


13.074 


IZ 


F 


Ft A T Oi^l 
[L/AijyuJ 




03 


29 


55.50 


+31 


15 


30.5 


15.120 


14.131 


IZ 


F 






03 


29 


52.65 


+31 


17 


22.9 


16.325 


14.970 


IZ 


F 






03 


29 


39.61 


+31 


17 


43.4 






JK 


F 






03 


28 


35.46 


+31 


21 


29.9 


16.052 


14.197 


IZ 


F 


Ft A T fli^l 
[L/AijyuJ 


4y 


03 


28 


48.45 


+31 


20 


28.4 


16.842 


14.283 


IZ 


F 


Ft A T rii^l 
[L/ALybJ 


100 


03 


29 


10.82 


+31 


16 


42.7 


15.652 


13.039 


JK 


F 


Ft A T fl^il 

[LALybJ 


235 


03 


29 


21.42 


+31 


15 


55.3 


15.396 


14.004 


IZ 


F 


[LAL96] 


305 


03 


29 


45.41 


+31 


16 


23.2 


15.118 


13.767 


IZ 


F 




03 


29 


50.23 


+31 


15 


47.9 


16.980 


15.106 


IZ 


F 






03 


30 


01.93 


+31 


10 


50.5 


14.374 


12.921 


IZ 


F 






03 


29 


54.78 


+31 


11 


41.7 


16.679 


15.253 


IZ 


F 






03 


29 


34.57 


+31 


11 


23.8 


16.750 


14.643 


IZ 


F 


[TAT Qfil 


OtJ-L 


03 


29 


36.00 


+31 


12 


49.6 


15.522 


14.615 


IZ 


F 


[LAL96] 


353 


03 


29 


26.11 


+31 


11 


36.9 


14.781 


12.874 


IZ 


F 


[LAL96] 


320 


03 


28 


01.19 


+31 


17 


36.5 


15.587 


14.136 


IZ 


F 




03 


27 


52.51 


+31 


19 


38.8 


14.973 


13.662 


IZ 


F 






03 


28 


22.90 


+31 


15 


21.7 


15.032 


13.587 


IZ 


F 


[LAL96] 


20 


03 


29 


08.71 


+31 


12 


01.9 






IZ 


F 






03 


29 


12.24 


+31 


12 


20.5 


16.588 


14.722 


IZ 


F 


[LAL96] 


254 


03 


29 


18.45 


+31 


11 


30.5 


16.413 


15.530 


IZ 


F 






03 


29 


31.00 


+31 


11 


20.1 


17.052 


15.265 


IZ 


F 






03 


29 


28.99 


+31 


10 


00.3 


13.365 


10.889 


IZ 


F 


[LAL96] 


329 


03 


29 


45.40 


+31 


10 


35.6 


14.828 


13.516 


IZ 


F 






03 


29 


46.48 


+31 


08 


43.6 


15.205 


14.096 


IZ 


F 






03 


29 


49.16 


+31 


09 


03.7 


16.143 


14.537 


IZ 


F 






03 


29 


28.95 


+31 


07 


40.9 


15.975 


15.199 


IZ 


F 


[LAL96] 


330 


03 


29 


24.73 


+31 


07 


26.8 






IZ 


F 






03 


29 


20.11 


+31 


08 


53.7 


15.923 


14.813 


IZ 


F 






03 


29 


09.23 


+31 


08 


55.4 


16.330 


14.912 


IZ 


F 






03 


28 


57.45 


+31 


09 


46.5 


16.623 


12.738 


IZ 


F 


[LAL96] 


160 


03 


28 


06.32 


+31 


10 


02.0 


15.731 


14.579 


IZ 


F 




03 


28 


22.07 


+31 


10 


42.9 


13.867 


11.859 


IZ 


F 


[LAL96] 


18 


03 


28 


39.01 


+31 


08 


07.6 


15.143 


13.389 


IZ 


F 


[LAL96] 


65 


03 


29 


03.60 


+31 


07 


11.9 


14.995 


13.351 


IZ 


F 


[LAL96] 


197 


03 


29 


05.23 


+31 


07 


10.6 


16.776 


14.766 


IZ 


F 





1 Photometry from 2MASS,if available 

^ Selected from IZ catalogue (IZ), JK plus Spitzer catalogue (JK) 
3 Source of s pectrum: M - MOIRCS, F - FMOS 

* [LAL96] - ILada etlHI (1996): MBO - Wilk ing et al.l I I2004I) : if no identifier is listed, the 



object does not have a known counterpart within 1" 



SONYC: NGC1333 



TABLE 7 

Objects excluded by spectroscopy by IScholz et al.1 ll2009Af) 



a(J2000) <5(J2000) J (mag)l K (mag)l Sel^ Spec^ Identifier-* 



03 


28 


41, 


.72 


+31 


11 


15, 


.1 










IZ 


M 




03 


28 


41, 


.97 


+31 


12 


17. 


.2 


15, 


.863 


13, 


.822 


IZ 


M 




03 


28 


46, 


.21 


+31 


12 


03, 


.4 


16. 


.933 


13, 


.526 


IZ 


M 


[LAL96] 90 


03 


28 


48, 


.99 


+31 


12 


45. 


.1 


17. 


.836 


14, 


.048 


IZ 


M 


[LAL96] 103 


03 


28 


52, 


.10 


+31 


16 


29. 


.3 


16. 


.071 


13, 


.738 


IZ 


M 


[LAL96] 123 


03 


28 


57, 


.25 


+31 


07 


26, 


.0 










IZ 


M 


[LAL96] 159 


03 


28 


58, 


.68 


+31 


09 


39, 


.2 


17. 


.052 


12, 


.945 


IZ 


M 


[LAL96] 169 


03 


29 


00, 


.70 


+31 


22 


00. 


.8 


16. 


.236 


11, 


.764 


IZ 


M 


[LAL96] 180 


03 


29 


17, 


.93 


+31 


14 


53. 


.5 


16. 


.625 


14, 


.039 


IZ 


M 


[LAL96] 287 


03 


29 


18, 


.66 


+31 


20 


17, 


.8 


17. 


.510 


14, 


.608 


IZ 


M 


MBO 109 


03 


29 


19, 


.86 


+31 


18 


47, 


.7 


17. 


.205 


13, 


.321 


IZ 


M 


[LAL96] 297 


03 


29 


28, 


.06 


+31 


18 


39. 


.0 


15. 


.307 


12, 


.846 


IZ 


M 


[LAL96] 327 


03 


29 


32, 


.20 


+31 


17 


07. 


.3 


15. 


.503 


13, 


.784 


IZ 


M 


[LAL96] 341 


03 


29 


37, 


.41 


+31 


17 


41, 


.6 


14. 


.886 


12, 


.761 


IZ 


M 


[LAL96] 359 


03 


29 


08, 


.17 


+31 


11 


54, 


.6 


17. 


.247 


14, 


.980 


IZ 


M 


[LAL961 217 



^ Photometry from 2MASS,if available 

^ Selected from IZ catalogue (IZ), JK plus Spitzer catalogue (JK) 
^ Source of spectrum: M - MO IRCS, F - FMOS 

^ [LAL96] - iLada et all ( fT996h ; MBO - iWilking et al.| | f200i) : if no identifier is listed, 



the object does not have a known counterpart within 1" 



16 



Scholz et al. 



Muzic, K., Scholz, A., Geers, V., Fissel, L., & Jayawardhana, R. 

2011, ApJ, 732, 86 
Natta, A., Testi, L., MuzeroUe, J., & et al. 2004, A&A, 424, 603 
Oppenheimer, B. R., Kulkarni, S. R., & Stauflfer, J. R. 2000, 

Protostars and Planets IV, 1313 
Padoan, P., Juvela, M., Kritsuk, A., & Norman, M. L. 2009, ApJ, 

707, L153 

Padoan, P., & Nordlund, A. 2004, ApJ, 617, 559 
Patience, J., King, R. R., de Rosa, R. J., & Marois, C. 2010, A&A, 
517, A76+ 

Preibisch, T., Stanke, T., & Zinnecker, H. 2003, A&A, 409, 147 
Quillen, A. C., Thorndike, S. L., Cunningham, A., & et al. 2005, 

ApJ, 632, 941 
Reipurth, B., & Clarke, C. 2001, AJ, 122, 432 
Scholz, A., Geers, V., Jayawardhana, R., & et al. 2009a, ApJ, 702, 

805 

Scholz, A., Jayawardhana, R., & Wood, K. 2006, ApJ, 645, 1498 
Scholz, A., Jayawardhana, R., Wood, K., & et al. 2007, ApJ, 660, 
1517 



Scholz, A., Xu, X., Jayawardhana, R., & et al. 2009b, MNRAS, 
398, 873 

Skemer, A. J., Close, L. M., Sziics, L., & et al. 2011, ApJ, 732, 107 
Slesnick, C. L., Hillenbrand, L. A., & Carpenter, J. M. 2004, ApJ, 
610, 1045 

Stamatellos, D., & Whitworth, A. P. 2009, MNRAS, 392, 413 
Umbreit, S., Burkert, A., Henning, T., Mikkola, S., & Spurzem, R. 

2005, ApJ, 623, 940 
Weights, D. J., Lucas, P. W., Roche, P. F., Pinfield, D. J., & 

Riddick, F. 2009, MNRAS, 392, 817 
Wilking, B. A., Meyer, M. R., Greene, T. P., Mikhail, A., & 

Carlson, G. 2004, AJ, 127, 1131 
Winston, E., Megeath, S. T., Wolk, S. J., & et al. 2009, AJ, 137, 

4777 

Zapatero Osorio, M. R., Bejar, V. J. S., Martm, E. L., & et al. 
2000, Science, 290, 103 



